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Abstract. 

Purpose: To improve the delivery efficiency of VMAT by extending the recently 
published VMAT treatment planning algorithm vmerge to automatically generate 
optimal partial-arc plans. 

Methods and materials: A high-quality initial plan is created by solving a convex 
multicriteria optimization problem using 180 equi-spaced beams. This initial plan is 
used to form a set of dose constraints, and a set of partial-arc plans is created by 
searching the space of all possible partial-arc plans that satisfy these constraints. For 
each partial-arc, an iterative flucnce map merging and sequencing algorithm (vmerge) 
is used to improve the delivery efficiency. Merging continues as long as the dose quality 
is maintained above a user-defined threshold. The final plan is selected as the partial- 
arc with the lowest treatment time. The complete algorithm is called pmerge. 
Results: Partial- arc plans are created using pmerge for a lung, liver and prostate 
case, with final treatment times of 127, 245 and 147 seconds. Treatment times using 
full arcs with vmerge are 211, 357 and 178 seconds. Dose quality is maintained across 
the initial, vmerge, and pmerge plans to within 5% of the mean doses to the critical 
organs-at-risk and with target coverage above 98%. Additionally, we find that the 
angular distribution of fluence in the initial plans is predictive of the start and end 
angles of the optimal partial-arc. 

Conclusions: The pmerge algorithm is an extension to vmerge that automatically 
finds the partial-arc plan that minimizes the treatment time. VMAT delivery efficiency 
can be improved by employing partial-arcs without compromising dose quality. Partial- 
arcs are most applicable to cases with non-centralized targets, where the time savings 
is greatest. 
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1. Introduction 

Volumetric modulated arc therapy (VMAT) is an extension to intensity modulated 
radiation therapy (IMRT) that offers the ability to deliver highly conformal dose 
distributions in a fraction of the time as IMRT (Yu & Tang 2011). In both techniques, 
the delivered radiation is modulated by a multi-leaf collimator (MLC) that conforms the 
dose to the targets. However, unlike IMRT where radiation is delivered at a fixed set 
of beam angles, in arc therapy the gantry delivers radiation continuously as it moves in 
an arc around the patient. In VMAT, one of the more widely used forms of arc therapy, 
the gantry speed, dose rate, and MLC leaf positions are allowed to vary as the gantry 
moves around the couch (Otto 2008). The additional degrees of freedom afforded by 
VMAT have been shown to maintain or improve the dose conformity when compared 
with IMRT (Teoh et al. 2011, Palma et al. 2008, Cozzi et al. 2008, Johnston et al. 2011). 
However, this also makes the treatment planning process much more complex. 

One of the most compelling reasons for using VMAT over IMRT is the ability to 
reduce the treatment time, often below 5 minutes per session (Bedford 2009, Clivio 
et al. 2009, Rao et al. 2010, Zhang et al. 2010). Reducing treatment time and the 
number of monitor units delivered reduces errors from patient motion (Hoogeman 
et al. 2008), decreases the scattered dose responsible for secondary malignancies (Hall 
& Wuu 2003, Ruben et al. 2008, Kry et al. 2005), and may have direct radiobiological 
benefits (Wang et al. 2003, Bese et al. 2007). In addition, it reduces the amount of time 
patients need to spend receiving their therapy, and may allow for the treatment of more 
patients. 

Approaches in reducing treatment time have focused on, for example, optimizing the 
delivery parameters (e.g. dose rate) given an initial plan (Rangaraj et al. 2010, Boylan 
et al. 2011), improving the delivery control system (Bertelsen et al. 2011) and developing 
leaf sequencing algorithms that improve delivery efficiency (Luan et al. 2008). In our 
own work on VMAT, we introduced vmerge, a method to facilitate the selection of 
the best balance between planning objectives and delivery efficiency, one of the main 
tradeoffs in VMAT treatment planning (Craft et al. 2012). Inherent to most approaches 
is the idea that in the tradeoff between dose quality and delivery efficiency, dose quality 
should be highly favored. 

Arc therapy can be delivered as a partial-arc, a full single arc, or as multiple arcs. 
For complex cases where a large amount of beam modulation is required, delivery in a 
single arc requires significant slowing of the gantry and reduction of the dose rate to 
achieve dose distributions similar to IMRT (Bortfeld & Webb 2009). Instead, radiation 
may be delivered over multiple arcs, which provides the necessary degrees of freedom 
but will tend to increase the total treatment time and required number of monitor 
units (Guckenberger et al. 2009). Alternatively, for cases with peripherally located 
targets, one might expect that appreciable beam modulation is only required over a 
critical angle range. In this case, the gantry might move at maximal speed and with small 
aperture sizes over the angles that do not provide for effective radiation. Eliminating 
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this portion of the arc from the delivery would be expected to reduce the treatment time 
with minimal effect on the dose distribution. However, to the best of our knowledge, no 
VMAT algorithm explicitly addresses this situation. In this work, we present a method 
called pmerge for automatically finding the optimal partial-arcs that eliminate the 
unneeded arc segments without significantly compromising dose quality. 

Like the beam angle optimization problem for fixed beam IMRT, the partial-arc 
optimization problem is nonconvex. Nonconvexity means that a local minimum is not 
guaranteed to be a global minimum. In order to guarantee a global minimum for 
nonconvex problems, a global search needs to be performed. Given that the nonconvex 
space of the partial-arc optimization problem consists of only two variables, the start 
and end gantry angles, this global search can be done by discretizing the search space 
and enumerating all possible start/end angle pairs, which is the strategy we adopt. This 
strategy is inherently parallelizeable, allowing for shorter computation times. 

To solve the VMAT problem for each partial-arc we use our previously reported 
merging and sequencing algorithm called vmerge (Craft et al. 2012). We call 
the complete process pmerge, which employs both merging and sequencing and an 
enumerable set of partial-arcs to reduce treatment time. Given an initial high quality 
plan as a dose "gold standard", pmerge will automatically determine the fluence 
maps, dose rate, gantry speed, and starting and ending gantry angles that minimize 
the treatment time while maintaining the dose distribution to within a user-defined 
quality threshold. 

2. Materials and methods 

We define the following terms: initial plan, partial- arc plan, planning goals, VMERGE 
plan and pmerge plan. The initial plan refers to a high-quality plan that serves as the 
dose standard used to constrain all partial-arcs. The initial plan is created without 
regards to delivery efficiency, providing a standard for understanding the tradeoff 
between dose quality and delivery efficiency Our strategy for creating the initial plan 
is described in Section 12.21 where we define a 180-beam IMRT problem, with beams 
placed every 2°, and solve it using multi-criteria optimization (MCO). 

Next, a set of partial-arc plans is created using a two-step process. First, we solve 
a set of restricted angle IMRT problems with various start and end gantry angles, 
with beams placed every 2°, by finding feasible solutions that minimize the sum-of- 
positive gradients (SPG), subject to matching the initial dose quality (Section [2.31) . We 
then use vmerge to iteratively merge neighboring fluence maps to create deliverable 
VMAT plans with improved delivery efficiency Merging is automatically halted when 
the merged plan violates a set of user-defined planning goals, which articulates how 
much the planner is willing to deviate from the initial plan in order to improve delivery 
efficiency We refer to the vmerge plan as the result of the merging algorithm on 
the 180-beam plan, and the pmerge plan as the merged partial-arc with the lowest 
treatment time. The full method is outlined in Figure [TJ 
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Figure 1. Illustration of our three step process. 1) An initial IMRT plan is created to 
generate a set of dose constraints for pmerge. 2a) pmerge uses these constraints to 
generate a set of partial-arc solutions. 2b) Each TV-beam partial-arc solution is merged 
and sequenced with VMERGE, subject to the user-defined planning goals. 3) The final 
plan is selected as the one with the lowest treatment time. 



2.1. Delivery parameters and leaf sequencing 

There are two classes of methods for solving VMAT problems: one-stage and two- 
stage (Yu &; Tang 2011). In the one stage method, known as direct aperture optimization, 
the paths of the MLC leaves are directly optimized, taking into account limitations of 
finite leaf speed. Instead, we use a two-stage planning approach to solve the VMAT 
problem for any size arc. In the first step, an iV-beam IMRT problem with beams 
spaced every 2° is solved to generate a set of N fluence maps {fl x , fi x ---, fi N x ), where 
each fluence map is a KxL matrix of L rows of MLC leaves and K columns. In the 
second step, a leaf sequencing algorithm determines the MLC leaf trajectories, gantry 
speed (cj), and dose rate (r) for each 2° arc portion, which we call a gantry sector. In 
leaf sequencing, the trajectories of the leaves are set so that the accumulated fluence 
over a given gantry sector matches the input fluence from the static beam fluence map. 
To accomplish this, we employ a widely used dynamic MLC (dMLC) sequencer which 
creates a "sliding window" that moves leaves unidirectionally across the field for each 
arc portion (Spirou & Chui 1994, Svensson et al. 1994). Thus, for each gantry sector, the 
leaves must move from left-to-right (or right-to- left), traversing the entire field of width 
W. In practice, the gantry speed, dose rate and maximum leaf speed are constrained 
by the limitations of the linac hardware and will have a significant effect on the total 
treatment time (Boylan et al. 2011). We select a set of hardware capabilities (Tabled]) 
that are representative of the values found on many current treatment machines. 

When applied to VMAT, the total treatment time using the dMLC algorithm is 
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Table 1. VMAT delivery parameters and constraints 



Delivery parameter 



Symbol Constraint 



Maximum leaf speed 
Maximum gantry speed 
Minimum gantry speed 



"max 2.5 cm/s 
w max 6 deg/s 
w min deg/s 



Maximum dose rate 
Minimum dose rate 



r max 600 MU/min 
r min 50 MU/min 



controlled by three terms: 1) the total distance the leaves must travel during delivery, 
2) the amount of beam modulation required to deliver the desired fluence map and 3) 
the maximum gantry speed. These three considerations can be combined to provide an 
equation to calculate the total treatment time for a plan with N gantry sectors of arc 
lengths A9f. 



where the (.) + operator is shorthand for max (0, .). The first term in (CQ) is the field 
width divided by the maximum leaf speed, which is the minimum time required for the 
leaves to move from the left side of the field to the right (or vice versa) if no fluence 
were delivered. As fluence is added, the leaf speed will be modulated to allow enough 
time for the dose to accumulate. This time is described in the second term as the 
sum-of-positive gradients (SPG), a measure of the amount of variation in a field (Craft 
et al. 2007), divided by the dose rate. Fields that are highly varied (have a high SPG 
term) will take longer to deliver with dMLC. Finally, because of limitations on the 
maximum gantry speed, the time required to deliver a gantry sector cannot be less than 
time required to move the gantry at maximum speed over the gantry sector. The total 
treatment time is the sum of the time required to deliver each of the N gantry sectors. 

Given a fixed set of hardware capabilities, the total treatment time can be improved 
by either decreasing the total number of gantry sectors N or by decreasing the SPG. 
Decreasing the number of gantry sectors motivates the partial-arc strategy described 
here, and the merging strategy vmerge - by using partial-arc solutions, we begin with 
a lower number of fluence maps than in a full arc, and this is decreased further by 
combing gantry sectors in the merging routine. Thus, for each gantry sector that is 
merged or eliminated with a partial-arc, the treatment time decreases by no less than 
W/v majK . Finally, to further reduce treatment time, we incorporate SPG minimization 
into our optimization step as described in Section 12.31 




+■ 
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2.2. Initial plan optimization 

We use a multicriteria optimization (MCO) approach to solve the initial 180-beam IMRT 
problem, as described in Craft et al. (2012). Briefly, in MCO, each structure is assigned 
an objective function, and a tradeoff surface of Pareto optimal plans is generated by 
the optimizer (Monz et al. 2008, Thieke et al. 2007). A Pareto optimal plan is one 
where no improvements could be made to any objective without worsening another 
objective (Cotrutz et al. 2001). The user navigates this tradeoff surface and selects the 
plan that best meets the treatment goals. 

The general multicriteria IMRT formulation is the following: 

optimize {g x (d) , g 2 (d), . . . g M (d)} 
subject to d = Df 
deC 

f > (2) 

where gi are the M objective functions describing the dose to a particular structure (e.g., 
mean dose to stomach), d is the vector of voxel doses, D is the dose-influence matrix, 
and / is a concatenation of all the fluence maps into a single beamlet fluence vector. 
The set C is a convex set of dose constraints (e.g. maximum dose to the target) that are 
designed to be met by all plans. We use the CERR 3.0 (Deasy et al. 2003) environment 
and its beamlet based quadratic infinite beam (QIB) dose computation to calculate the 
dose influence matrix, with a computation time of around 10 minutes per case. Given a 
dose matrix, the dose distribution can be calculated from the fluence maps in seconds. 

The constraints on the planning target volumes (PTV) and the OARs for each site 
are listed in Table [2j In our MCO formulation, we had the optimizer minimize the 
mean doses to the OARs. Additionally, we include as an OAR unclassified tissue (u.t.), 
which encompasses all voxels that do not belong to another structure. Optimizations are 
performed using the solver described by Chen et al. (2010). We then choose a plan for 
each treatment site by exploring the Pareto surface and selecting a plan that balanced 
the dose the OARs. 

Table 2. Target coverage constraints and OAR list used in the MCO optimization 



Treatment site PTV Min; Max (Gy) OARs 



Lung 


50.0: 


55.0 


left lung, right lung, esophagus, 
heart, spinal cord, u.t. 


Liver 


50.0: 


55.0 


left lung, right lung, liver, spinal cord, 
stomach, left kidney, right kidney, u.t. 


Prostate 


79; 5 


S5.3 


left femoral head, right femoral head, 
anterior rectum, bladder, u.t. 
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2.3. Generating the partial-arc plans 

The set of partial-arc solutions are created by iteratively solving the 2° grid IMRT 
problem for a range of partial-arcs. In creating the plans, we respect the constraint 
that the gantry cannot rotate through 0°, a line from the couch to the floor. Plans are 
created with arc lengths every 40° (360°, 320°, ...) until there are no feasible solutions. 
For each arc length, plans are created with arcs centered at intervals of 40°. Thus, there 
is one 360° plan, three 320° degree plans (beginning at — 40°, 0°, 40°), five 280° plans, 
etc. 

Each partial-arc solution is created to maintain the same dose quality as the initial 
180-beam solution. This is done with the same solver used to solve the original 180- 
beam MCO formulation, except run in feasibility mode, where the solver returns the 
first solution it finds which satisfies all of the hard constraints. Because treatment time 
scales with the complexity of the solution, this solution is then smoothed using total 
SPG as the objective function, while maintaining dosimetric feasibility: 



SPG minimization subject to dose matching 



minimize 



N 

y max 

i=1 rows I 




dx 



subject to d — Df 

di> B L) Vi G target 
di < Bu, Vi 

meanf\d) < mean^d), V OAR j 
/>0 



(3) 



where B L and B v is the target lower bound, By is the upper bound on the dose, and 



mean^(d) and mean^d) are the mean doses to the jth OAR for the partial- arc and 
initial plans, respectively. This formulation ensures that for each partial-arc solution 
the target coverage is no less than, and the maximum dose and the mean doses to each 
OAR are no greater than those of the initial plan. 



.(■0 



2.4- Merging and final plan selection 

After the SPG minimization step, we have a set of IMRT solutions for various partial- 
arcs. To improve the delivery efficiency, we use a previously described fluence map 
merging strategy called VMERGE to combine gantry sectors (Craft et al. 2012). Briefly, 
VMERGE employs a greedy merging strategy to iteratively add similar neighboring 
fluence maps - for each iteration, the two neighboring fluence maps that are most 
similar are added together into a single fluence map. This reduces the number of gantry 
sectors and improves the delivery efficiency, at the expense of reducing the amount of 
fluence modulation and the angular resolution. 
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pmerge automates vmerge by utilizing a quality cutoff, which limits how much 
the dose can change, and automatically stops merging when the quality drops below the 
planning goals. Here, the planner is free to chose how much the dose quality is allowed 
to change in order to improve delivery efficiency We use the initial high-quality plan 
as a guide for how to set these goals. For the OARs that receive the most radiation, we 
limit the partial-arc mean dose to less than 105% of the initial mean dose for that OAR. 
Additionally, we ensure that the target coverage stays above 98%, and the maximum 
dose is less than 102% of the initial maximum dose. Finally, we illustrate the use of 
site-specific planning goals, including limiting the maximum dose to the femoral heads 
in the prostate case and adding additional constraints on the left lung dose in the lung 
case. The full sets of planning goals are provided in Tables 3-5. 

3. Results 

3.1. Lung 



Table 3. Lung: Planning goals, dose quality and treatment times 





Goal 


Initial (I) 


PMERGE (P) 


u.t. dose 


mean[^ < 1.10 x mean^ 


2.8 Gy 


3.0 Gy 


Right lung 


mean^2 < 1-05 x mean^ 


0.6 Gy 


0.5 Gy 


Left lung 


t\/40( p ) < 1.05 x 


9.4% 


9.5% 




V50( p ) < 1.05 x 


3.0% 


3.0% 


Target coverage 


V50( p ) > 0.98 x V5QW 


100% 


98.2% 


Maximum dose 


max( p ) < 1.02 x max' 7 ' 


55.0 Gy 


55.8 Gy 


(Start, end) angles 




(0°, 360°) 


(60°, 260°) 


Treatment time* 




211 s 


127 s 



\Vx denotes the volume receiving at least x Gy 
^Treatment times are after merging 



We navigate to an initial 180-beam plan with mean doses and target coverage as 
shown in Table [3J The treatment time of the initial plan before merging is 771 s, and 
is reduced to 211 s with vmerge. The initial plan and vmerge plans are shown with 
the outer two rings of Fig [2(c). The VMERGE plan is notable for a large arc portion of 
116° on the contralateral side of the PTV where the gantry speed is at its maximum. 
As shown in Figure [2(d), there is little fluence delivered over this section in the initial 
180-beam plan. 

Partial-arc solutions were attempted for arcs with lengths between 80° and 360°, at 
intervals of 40°. 50 had feasible solutions that met the dose constraints from Equation [31 
and 26 of these had lower treatment times than the vmerge plan. The final pmerge 
plan has an arc length of 200° and treatment time of 127 s, and is shown in the inner 
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two rings of Fig 2c, both before and after merging. The total difference in treatment 
times between the vmerge (211 s) and pmerge (127 s) plans is 84 s, a 40% reduction. 
Optimization, merging and sequencing took about 40 minutes to complete all 64 plans, 
running in parallel on 8 processors. Computation time was similar across the lung, liver 
and prostate cases. 

3.2. Liver 

Table 4. Liver: Planning goals, dose quality and treatment times 



Goal Initial (I) pmerge (P) 



u.t. 


mean|^ < 1.10 x mean^ 


5.7 Gy 


5.6 Gy 


Right lung 


mean^2 < 1-05 x mean^ 


9.1 Gy 


9.2 Gy 


Left lung 


no constraint (mean^) 


1.5 


1.5 Gy 


Liver 


mean^ < 1.05 x mean^ 


23.3 Gy 


23.3 Gy 


Right kidney 


mean^- < 1.05 x mean^ 


5.8 Gy 


4.7 Gy 


Left kidney 


( p) 

no constraint (mean^) 


0.7 Gy 


0.4 Gy 


Target coverage 


V50^ > 0.98 x 


100% 


98.0% 


Maximum dose 


max( p ) < 1.02 x max'^ 


55.0 Gy 


55.6 Gy 



(Start, end) angles (0°, 360°) (100°, 340°) 



Treatment time^ 357 s 245 s 

f Treatment times are after merging 

We navigate to an original 180-beam solution with mean OAR doses, target 
coverage and maximum dose as shown in Table HI The treatment time before merging 
for the full arc solution is 1214 s and is reduced to 357 s with vmerge (Figure EJ^c), 
outer rings). Partial-arc solutions were attempted for arcs lengths between 80° and 
360°, at intervals of 40°, and 33 had feasible solutions. There were no solutions with arc 
lengths less than 160°. Of the feasible partial-arcs, 19 had lower treatment times than 
the vmerge plan. 

The selected partial-arc has an arc length of 240° (Figure EJ^c), inner rings) and 
treatment time of 245 s. The total difference in treatment times between the vmerge 
plan (357 s) and pmerge plan (245 s) is 112 s, a 31% reduction. In this case the initial 
180-beam solution is significant for two peaks of fluence, one at approximately 150° and 
one at 330°, as shown in Figure [3]^d). The optimal partial-arc includes both of these 
peaks, and is able to eliminate the first 100° where the amount of required fluence is 
low. 
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Table 5. Prostate: Planning goals, dose quality and treatment times 



Goal Initial (I) pmerge (P) 



u.t. 


rpi m — 

mean^ < 1.10 x meanly 


15.5 Gy 


16.4 Gy 


Anterior rectum 


(P) 1 ^i- (I) 
mean^ < 1.05 x mean^jj 


41.4 Gy 


41.8 Gy 


Bladder 


mean^ < 1.05 x mean^ 


39.6 Gy 


39.8 Gy 


Left femoral head 


max^ < 1.10 x ma.x'lpfj 


33.6 Gy 


35.5 Gy 


Right femoral head 


max^ H < 1.10 x m&Xftp H 


37.0 Gy 


38.3 Gy 


Target coverage 


V79^ > 0.98 x V79& 


100% 


98.2% 


Maximum dose 


max' ? ' < 1.02 x max'^ 


85.3 Gy 


86.7 Gy 


(Start, end) angles 




(0°, 360°) 


(40°, 320°) 


Treatment time* 




178 s 


147 s 



f Treatment times are after merging 



3.3. Prostate 

We navigated to an original 180-beam solution with mean OAR doses, target coverage 
and maximum dose as shown in Table [51 The treatment time before merging for the full 
arc solution is 860 s and is reduced to 178 s with vmerge (Figure H^c), outer two rings). 
Unlike the lung and liver cases, the vmerge plan has regularly spaced arc portions of 
similar size, and the fluence is delivered fairly uniformly around the target, as shown in 
Figure |4(d). This is reflected in the optimal partial-arc solution, which is highly similar 
to the initial plan through nearly the full arc. 

Partial-arc solutions were created for arcs lengths between 80° and 360°, at intervals 
of 40°, and 45 had feasible solutions. Of the 45, only 3 had treatment times lower 
than the vmerge plan. There were no feasible partial- arcs at 80°. For plans with 
arc lengths of less than 280°, almost no merging was possible before the plan quality 
degraded beyond the stopping criteria. The selected partial-arc has an arc length of 
300° (Figure H](c), inner rings) and treatment time of 147 s. The total difference in 
treatment times between the vmerge (178 s) plan and pmerge plan (147 s) is 31 s, a 
17% reduction. 

4. Discussion and conclusions 

We present a method for automatically creating efficient VMAT partial-arc plans from 
a high quality reference plan. By sampling a set of partial-arcs, we reduced treatment 
times for a lung, liver and prostate case. The final pmerge plans reflect the underlying 
geometry of the problem. For the lung and liver cases with non-centralized PTVs, a 
signifiant portion of the full-arc can be eliminated with little effect on the dose quality. 
This allowed for the treatment time to be reduced by 40% and 31% respectively. The 
reduction in treatment time for the prostate case was more modest at 17%. As might be 
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expected for a centrally located target, the fTuence was distributed more evenly around 
the arc, and most of the partial- arcs were not able to improve the delivery efficiency 
while maintaining the dose quality to within the specified boundaries. 

In pmerge, the dose quality of the final partial-arc is designed to be highly similar 
to the initial plan. However, the pmerge plan is limited only by the planning goals 
specified after the initial plan creation. Thus, the partial-arc plan could deviate from 
the initial plan in ways that are not specified by this quality threshold. Any number 
or type of constraints could easily be imposed on the final solution to ensure that the 
initial plan and pmerge plan were sufficiently similar. In practice, we did not find that 
it was necessary to impose many constraints in order to match the solutions. Although 
we were free to use any number of dose-volume points or EUD functions, our planning 
goals were restricted to only one or two constraints per OAR, usually a limit on the 
mean or maximum dose. For all three cases, these were sufficient to closely match the 
doses across the DVH, even for OARs which were not explicitly specified. 

In VMAT, as in IMRT, there are many different ways to deliver the fluence that 
will give highly similar dose distributions. Because the treatment time is dependent 
on the geometry of these fluence maps, the method by which the initial fluence maps 
are created is an important determinant of treatment time. Although pmerge is not 
confined to any single optimization strategy, we found that using SPG smoothing was 
critically important to achieving highly efficient plans. Indeed, fluence map smoothing is 
important in IMRT (e.g. (Matuszak et al. 2008)), where it is know that the complexity 
of the fluence maps - and hence the number of monitor units - can be greatly reduced 
without significantly impacting the dose distribution. 

Although the partial-arc global search can be parallelized, we recognize that the 
high computation time is a limit to our method. The observation that the optimizer 
distributes fluence very similarly between the full-arc and partial-arc plans suggests a 
way of approximating the appropriate partial-arc without having to enumerate all of the 
possibilities. One strategy would be to plot the arrangement of fluence around the arc 
and have the planner select the partial-arc that captures the critical angles. For the lung 
and liver cases, we find that there were many different partial-arcs that greatly improved 
the treatment time, and thus it is not essential to find the absolute global minimum 
when plans near the minimum are almost equivalent (see Figures [2(a) and 131(a)). More 
importantly, we conclude that partial-arcs that cover certain critical angles are useful 
for reducing treatment time in VMAT, and pmerge is a useful first-pass method at 
finding the optimal plan. 
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Figure 2. Lung: a) Final treatment times for each of the merged partial-arc plans 
(blue squares) with treatment times near or below the vmerge plan (red circle), b) 
DVH for the initial 180-beam plan (solid) and the PMERGE plan (dashed), c) Initial 
and merged gantry sectors for the pmerge plan (inner rings) and 180-beam plan (outer 
rings) over the pmerge dose wash, d) Distribution of flucnce for the initial solution 
and the optimal partial-arc before merging (100-beam). 
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Figure 3. Liver: a) Final treatment times for each of the merged partial-arc plans 
(blue squares) with treatment times near or below the vmerge plan (red circle), b) 
DVH for the initial 180-beam plan (solid) and the PMERGE plan (dashed), c) Initial 
and merged gantry sectors for the PMERGE plan (inner rings) and 180-beam plan (outer 
rings) over the PMERGE dose wash, d) Distribution of fluence for the initial solution 
and the optimal partial-arc before merging (120-beam). 
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Figure 4. Prostate: a) Final treatment times for each of the merged partial-arc 
plans (blue squares) with treatment times near or below the vmerge plan (red circle) . 
b) DVH for the initial 180-beam plan (solid) and the pmerge plan (dashed), c) Initial 
and merged gantry sectors for the pmerge plan (inner rings) and 180-beam plan (outer 
rings) over the PMERGE dose wash, d) Distribution of fluence for the initial solution 
and the optimal partial-arc before merging (140-beam). 



